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We present results at — 6.0 and 6.2 for the O(a) improvement and renormalization constants for bilinear 
operators using axial and vector Ward identities. We discuss the extraction of the mass dependence of the 
renormalization constants and the coefficients of the equation of motion operators. 
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1. INTRODUCTION 

In quenched Lattice QCD, axial and vector 
Ward identities can be used to determine, at 
0(a), all the scale independent renormalization 
constants for bilinear currents [Za, Zy, and 
Zs/Zp), the improvement constants {ca, cy, and 
ct), the quark mass dependence of all five Zq, 
and the coefficients of the equation of motion op- 
erators p]-|(| ■ Here we summarize results at (3 — 
6.0 and 6.2 and discuss the highlights of our cal- 
culations. Details and notation are given in B. 

We start with the generalized axial Ward iden- 
tity involving operators improved on and off-shell 

(6S^ o^%(y) J^(0)) = (60^ ff (y) J^(0)) 

where SO is the result of the axial variation of O 

), and SS 



P, and T H 



is the variation in the action. 

At 0{a) there exists only one dimension 4 off- 
shell operator (which vanishes by the equations of 
motion) for each bilinear that has the appropriate 
symmetries M. Consequently, we define 



O 



m 

R,off 

/n(y) 
U I,off 



y(ij) fl('i) 

W I,off ' 



E, 



(ij) 



0^' (1/4) ac' Eg j) , 



(1) 

(2) 
(3) 



Here (ij) (with i ^ j) specifies the flavor, and 
V$ipj = C0 + m^ijjj + 0(a 2 ) is the full O(a) im- 
proved Dirac operator for quark flavor j. This 
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ensures that the equation-of-motion operator Eq 
gives rise only to contact terms, and does not 
change the overall normalization Zq. The O(a) 

(ij) 

on-shell improved renormalized operators O r 
are 



O 



(ij) 



7° 



(1 



= Z^l + boam^O} 



(ij) 



(Vi), 
Pi 

The Z9 



= A, 

= V a 



= P 



ft V 

P, 



- ac A d^P , 
acvdvT^v , 
f ac T (d^V v - 
Si = S j 



duV n 



(4) 

(5) 
(6) 
(7) 
(8) 
(9) 



^ are renormalization constants in the chi- 
ral limit, rriij = (m, + rrij)/2 is the average bare 
quark mass, ami = 1/2/tj — 1/2k c , k c is the value 
of the hopping parameter in the chiral limit, and 
my is the quark mass defined by the axial Ward 
identity (AWI) in Eq. [l^. Note that m and m are 
identical in a discretized theory with chiral sym- 
metry, like staggered fermions. With these defi- 
nitions, bo = 1, cq = 0, c' — 1 at tree level Q. 

Since the equation-of-motion operators con- 
tribute only contact terms, Eq. (|l]) can be rewrit- 
ten in terms of just on-shell improved operators: 

(Ld^xSSiOf^iy^ff) J( 31 )(0)) 



{SOf\y 4l y) J(3D(0)> 



7 (13) 



,(12) 7 (23) 



—mi2 



0(a 2 ) (10) 



where 

SS I (x) = 2rh 12 P i I 12) (x) - d^Aj^ix) 
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Table 1 

Simulation parameters, statistics, and the time interval in X4 defining the volume V over which the 
chiral rotation is performed in the AWL The lattice spacing is fixed using tq = 0.5 fermi, and is thus 
independent of the fermion action. The source J is placed at t = 0. 



Label 


P 


csw 


a- L (GeV) 


Volume 


L (fm) 


Confs. 


X4 


60TI 


6.0 


1.4755 


2.12 


16 3 x 48 


1.5 


83 


4-18 


60NPf 
60NPb 


6.0 


1.769 


2.12 


16 3 x 48 


1.5 


125 
112 


4-18 
27-44 


62NP 


6.2 


1.614 


2.91 


24 3 x 64 


1.65 


70 
70 


6-25 
39-58 



Our calculation is limited to the case fh\ — rri2 
(this simplification was used due to limited com- 
puter resources), in which case the r.h.s. of Eq. |l^ 
reduces to 



+ {bso - be 
(bso ~ bo) 



arris 



b A 



arhi (11) 



7° 
^so 

A 

7° 
^SO 

A ^O 

where rhi = rhij \ mj=mi . Using Eqs. [w| and [H], all 
the bo (except 6t which requires m\ ^ m-i), Co, 
c'q , and the scale independent normalization con- 
stants are determined by making suitable choices 
for J, O, and y in Eq. |l0| and studying it as a 
function of rh\ and 7713 (Eq. |Tl| ). 



2. RESULTS 

The lattice parameters used in our calculation 
are given in Tab. [I] and & Our final results, which 
supercede those in Ref. |5J| , are given in Tab. § 

In many cases a given on-shell improvement 
and normalization constants can be determined 
in a number of ways as discussed in || . Results 
in Tab. || are based on the AWI with the best sig- 
nal and smallest error. Table ||| also includes re- 
sults by the ALPHA collaboration Jl|-§] and the 
one-loop tadpole improved perturbative results. 
To simplify comparison with previous results, we 
quote both by, b A and by, b A - 

One of the goals of our calculation is to quantify 
the residual 0(a 2 ) errors and to understand the 
shortcomings of 1-loop perturbation theory. For 
0(a 2 ) errors we use two estimates: (i) the differ- 
ence between our results and those by the AL- 
PHA collaboration |l|-||], and (ii) the difference 



between using 2-point and 3-point discretization 
of the derivatives p] in the extraction of c A from 



£ 2 <^ + ac^P] W (x, t)J&) (0)) 



E 5 (^ (ii) (^*)J (ii) (o)> 



= 2m 



i] 1 



(12) 



and the subsequent effect of the difference in c A 
on other constants. This latter variation is quoted 
as the second error in Tab. |. 

These differences are compared to the ex- 
pected size of the residual discretization errors: 
(qAqcd) ~ 0.15 and 0.1 for the improvement 
constants and (oAqcd) 2 ~ 0.02 and 0.01 for the 
normalization constants at /3 = 6.0 and 6.2 re- 
spectively. 

A comparison, at /3 = 6.0, between simulation 
at csw = 1-4755 (tadpole improved theory) and 
csw = L769 (non-perturbatively 0(a) improved 
theory) shows that all the constants are sensitive 
to the choice of csw- It is therefore important to 
use csw determined non-perturbatively. 

The most significant comparison is between our 
results and those of the ALPHA collaboration. 
The only results which do not agree within 2-er 
statistical errors are those for Zy, c A and cy at 
(3 = 6, and for Z v at (3 = 6.2. The differences for 
Z v arc of size 0.01 and 0.005 at p = 6 and 6.2 
respectively, and are thus consistent with the ex- 
pected differences of 0(a 2 ). The differences for c A 
and cy are also consistent with the size expected 
of 0(a) differences, but are more notable because 
they correspond to very large fractional differ- 
ences (e.g. our c A at P = 6 has less than half the 
magnitude of that found by the ALPHA collabo- 
ration) . What we learn is that (i) co , which van- 
ish at tree level and are numerically small, depend 
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Table 2 

Values of k used in the three simulations, and the corresponding values of aM^ and the quark mass am 
extracted, m is defined by the AWI in Eq. [l^. n c is the zero of m obtained from quadratic fits in 1/k. 
The non-zero value of aM n at k c is indicative of the inadequacy of quadratic fits, a 2 M 2 as a function of 
l/2n, used to extract it, and discretization errors. Of these, the first is the dominant cause and points 
to the need for including quenched chiral logs in the fits [||. 







60TI 






60NP 






62NP 




Label 


K 


am 


aM n 


K 


am 


aMn 


K 


aih 


aM r 


Kl 


0.11900 


0.443(8) 


1.530(1) 


0.1300 


0.144(1) 


0.711(2) 


0.1310 


0.1345(6) 


0.609(1) 


K2 


0.13524 


0.105(1) 


0.571(2) 


0.1310 


0.118(1) 


0.630(2) 


0.1321 


0.1054(4) 


0.522(1) 


K3 


0.13606 


0.084(1) 


0.504(2) 


0.1320 


0.092(1) 


0.544(2) 


0.1333 


0.0727(3) 


0.418(1) 


K4 


0.13688 


0.063(1) 


0.431(2) 


0.1326 


0.075(1) 


0.488(2) 


0.1339 


0.0560(2) 


0.360(2) 


K 5 


0.13770 


0.042(1) 


0.348(3) 


0.1333 


0.056(1) 


0.416(2) 


0.1344 


0.0419(2) 


0.307(2) 


kq 


0.13851 


0.020(1) 


0.244(4) 


0.1342 


0.032(1) 


0.308(3) 


0.1348 


0.0306(2) 


0.261(2) 




0.13878 


0.013(1) 


0.195(8) 


0.1345 


0.025(4) 


0.262(12) 


0.1350 


0.0248(1) 


0.235(2) 


K c 


0.13926(2) 





0.082(15) 


0.13532(3) 





0.083(20) 


0.135861(5) 





0.066(10) 



substantially, at (3 = 6, on the method/definition 
used to extract them; (ii) the variation between 2- 
pt and 3-pt derivatives significantly smaller than 
the difference between our results and those of the 
ALPHA collaboration; and (iii) these differences 
in cy, and even more so in ca, are substantially 
reduced at (3 = 6.2. The change appears too rapid 
to be an O(a) effect. 

Both cy and ct are obtained as a small differ- 
ence between two large terms. Nevertheless, we 
are able to design Ward identities that yield these 
quantities with reasonable precision. In particu- 
lar, the significant improvement we obtain in de- 
termining cy using methods described in || re- 
duces the error in Z A , Z p jZ Q s , ct and c' A as the 
uncertainty in cy feeds into these quantities. 

When comparing against perturbative esti- 
mates the yardstick we use for the missing higher 
order terms is ~ a 2 s w 0.02 and 0.016, respec- 
tively. We find that tadpole-improved 1-loop per- 
turbation theory underestimates the deviations 
of renormalization and improvement constants 
from their tree level values. In all but one case, 
however, these discrepancies can be understood 
as a combination of a 2-loop correction of size 
(1 - 2) x oi 2 s [for Z v , Z% and c A ], higher or- 
der discretization errors of size (1 — 2) x oAqcd 
[for cy, ct and by], and statistical errors [for b A , 
bp, and bs]. The only exception is Z P jZ Q s , for 
which a very large higher order perturbative con- 



Table 4 



Results for off-shell mixing coefficients. 





60TI 


60NPf 


60NPb 


62NP 




+3.72(73) 


+2.38(50) 


+3.00(37) 


+1.72(16) 


c'a 


+3.28(94) 


+1.99(56) 


+2.45(46) 


+1.53(20) 




-0.98(76) 


+0.44(49) 


-0.33(29) 


+0.91(12) 


c's 


+3.00(73) 


+2.00(48) 


+2.72(33) 


+1.49(14) 




+3.24(75) 


+1.96(49) 


+2.60(38) 


+1.51(15) 



tribution of size 4 x a 2 s is needed to reconcile our 
non-perturbative results with 1-loop perturbation 
theory. 

In Tab. ||, we present, first results for the 
equation-of-motion improvement constants c' x . 
The combination c' P + c! is extracted by study- 
ing the dependence of Eq. |10| o n rh\ once the 
other constants defined in Eq. |ll|have been deter- 
mined. The errors in the determination of the c' 
are dominated by two quantities: (i) The uncer- 
tainty in ca feeds into the extraction of c' A , and 
(ii) the correlation function from which c' P + c' P 
is determined has a poor signal (the intermediate 
state is a scalar for J = S, O = P and SO = S 
in Eq. |To| ) . The uncertainty in c' p then feeds into 
c' v , c' s , and c' T . Overall, we find a very significant 
improvement in the quality of the results with in- 
creasing f3, i.e., between (3 — 6.0 and (3 = 6.2. 

Finally, we comment on results presented in 
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Table 3 

Final results for improvement and renormalization constants. The first error is statistical, and the 
second, where present, corresponds to the difference between using 2-point and 3-point discretization of 
the derivative used in the extraction of ca- The ALPHA collaboration results are from 0-0]. For the 
tadpole improved perturbative results (labelled P.Th.) see appendix in Ref. |(| and references there. 





= 6.0 


(3 = 6.2 




LANL 


LANL 


ALPHA 


P. Th. 


LANL 


ALPHA 


P. Th. 




1.4755 


1.769 


1.769 


1.521 


1.614 


1.614 


1.481 


V 


+0.747(1) 


+0.770(1) 


0.7809(6) 


+0.810 


+0.7874(4) 


+0.7922(4)(9) 


+0.821 


Z°A 


+0.791(7)(4) 


+0.807(2)(8) 


0.7906(94) 


+0.829 


+0.818(2)(5) 


+0.807(8) (2) 


+0.839 


Z P /Z% 


+0.811(9)(5) 


+0.842(5)(1) 


N.A. 


+0.956 


+0.884(3)(1) 


N.A. 


+0.959 


CA 


-0.022(6)(1) 


-0.037(4)(8) 


-0.083(5) 


-0.013 


-0.032(3)(6) 


-0.038(4) 


-0.012 


cv 


-0.25(5)(3) 


-0.107(17)(4) 


-0.32(7) 


-0.028 


-0.09(2)(1) 


-0.21(7) 


-0.026 


C T 


+0.09(2)(1) 


+0.06(1)(3) 


N.A. 


+0.020 


+0.051(7)(17) 


N.A. 


+0.019 


bv 


+ 1.44(3)(2) 


+1.43(1)(4) 


N.A. 


+ 1.106 


+1.30(1)(1) 


N.A. 


+1.099 




+ 1.53(2) 


+1.52(1) 


+ 1.54(2) 


+ 1.274 


+1.42(1) 


+ 1.41(2) 


+1.255 


b~A - b v 


-0.51(9)(4) 


-0.26(3)(4) 


N.A. 


-0.002 


-0.11(3)(4) 


N.A. 


-0.002 


bA - b v 


-0.49(9)(4) 


-0.24(3)(4) 


N.A. 


-0.002 


-0.11(3)(4) 


N.A. 


-0.002 


bp - b s 


-0.07(9)(2) 


-0.06(4)(3) 


N.A. 


-0.066 


-0.09(2)(1) 


N.A. 


-0.062 


bp - b A 


-0.126(58)(1) 


-0.07(4)(5) 


N.A. 


+0.002 


-0.09(3)(3) 


N.A. 


+0.001 


b A 


+0.92(10)(6) 


+1.17(4)(8) 


N.A. 


+ 1.104 


+1.19(3)(5) 


N.A. 


+1.097 


b A 


+ 1.05(9)(4) 


+1.28(3)(4) 


N.A. 


+ 1.271 


+1.32(3)(4) 


N.A. 


+1.252 


bp 


+0.80(11)(6) 


+1.10(5)(13) 


N.A. 


+ 1.105 


+1.11(4)(7) 


N.A. 


+1.099 


bs 


+0.87(14)(4) 


+1.16(6)(11) 


N.A. 


+ 1.172 


+1.19(4)(6) 


N.A. 


+1.161 



two recent papers. Using the Schrodinger func- 
tional, Ref. ||] calculates bA — bp, bs and 
Z P /{Z° S Z A ) for a range of (3 > 6. Their most 
striking result is that different discretizations of 
derivatives lead to very different results for b A — 
bp. For example, at (3 — 6, this quantity varies 
roughly from 0.17 to —0.17. While our number 
lies within this range, our estimate of 0(a) un- 
certainties is clearly a substantial underestimate. 

Reference |)| has determined ca using the same 
method and similar lattice parameters as here 
but with significantly more configurations. They 
study, at one k (~ k$ at both (3 = 6.0 and 6.2), 
the effect of using derivatives that are tree-level 
improved through 0(a 2 ) (our 3-pt), 0(a A ) and 
0(a 6 ). They find a larger dependence than what 
we get between 2-pt and 3-pt discretizations at 
K5. The 0(a 2 ) errors in the two calculations are, 
however, different due to the choice of source and 
the fit range in time. Also, we find that after chi- 
ral extrapolation these discretization effects are 



significantly reduced. Nevertheless, once again 
the large variation should serve as a warning that 
the O(a) errors in cq can be substantial. 

REFERENCES 



M. Liischer et al., Nuc. Phy. B478 (1996) 365. 
M. Liischer et at, Nuc. Phy. B491 (1997) 323. 
M. Liischer et al, Nuc. Phy. B491 (1997) 344. 
M. Bochicchio et a/., Nuc. Phy. B262 (1985) 
331. 

T. Bhattacharya et al, Phys. Lett. B461 
(1999) 79 

T. Bhattacharya et al, |hep-lat/0009038j . 

G. Martinelli et al, Phys. Lett. B411 (1997) 

141 



M. Guagnelli et al, |hep-lat/0009021|. 

S. Collins and C. Davies, |hep-lat/0010045 



